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, The pulse control of decoherence in a qubit interacting with a quantum environment is 



studied with focus on a general case where decoherence is induced by both pure dephasing 



I and population decay. To observe how the decoherence is suppressed by periodic tt pulses, 

Q I we present a simple method to calculate the time evolution of a qubit under arbitrary pulse 

I sequences consisting of bit- flips and/or phase-flips. We examine the effectiveness of the two 

typical sequences: bb sequence consisting of only bit-flips, and bp sequence consisting of both 

' bit- and phase-flips. It is shown that the effectiveness of the pulse sequences depends on a 

, relative strength of the two decoherence processes especially when a pulse interval is slightly 
00 ' 

o, 

OO bp sequence is always more effective than, or at least as effective as, the bb sequence. 

o ■ 

*n ' KEYWORDS: decoherence, pure dephasing, population decay, qubit 
O 

:^ 

C ' 1. Introduction 

^ . To employ a two-level quantum system as a qubit for quantum computations, a quantum 



superposed state prepared in it must be well preserved for a long time. However, since any 
quantum systems interact with their surrounding environment, a superposed state is degraded 



^ I by getting entangled with an environment. The degradation of a quantum superposed state 

is called decoherence. To perform reliable quantum computations, we thus need to suppress 
decoherence caused by an environment. Several schemes to suppress decoherence have been 
proposed so far.^"^ Among them, we focus on the decoherence suppression scheme which uses 
a sequence of short vr pulses^' ^'^"^^ since it is relatively easy to implement in experimental 
situations and can be readily applied to multi-qubit systems. We call this scheme the pulse 
control of decoherence. The pulse control has been proposed by Ban,^ and by Viola and Lloyd. ^ 
They considered a qubit coupled with an environment assuming that qubit-environment inter- 
actions induce decoherence without population decay (i.e., pure dephasing). They calculated 
the time evolution of a qubit under periodic vr pulses, and showed that the pure dephasing 
can be suppressed if the pulse interval is much shorter than the correlation time for qubit- 
environment interactions. It has been proposed that we can even suppress the decoherence 
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with population decay by applying vr pulses.^' "^'^ 

Let us consider a qubit linearly interacting with a boson environment. The Hamiltonian 
is Ho = Hs+'Hb+ Hsb with 

Hs = /li^oy^lB, (1) 
Hb = ls^J2^kblbk. (2) 

k 

Here, T^s and Hb describe the qubit and the environment, respectively, and Hsb describes 
their mutual interactions. We have used the pseudo-spin representation in expressing the 
qubit. Ban, and Viola and Lloyd, considered the case of T^sB = 'Hz with 

Hz = ^fi [plz^^z ®h\ + QkzCFz <8) bk^ , (3) 

k 

which induces the pure dephasing.^'^ They pointed out that the dephasing can be suppressed 
by the sequence of periodic tt pulses about the x-axis in the pseudo-spin space. They confirmed 
that the dephasing is notably suppressed if the interval At between tt pulses is much shorter 
than the correlation time for qubit-environment interactions. Uchiyama and Aihara^^ extended 
the argument to a nonlinear interaction case by using a diagrammatic perturbation method. 
However, we should note that the decoherence arises from not only pure dephasing but also 
population decay. A simple but standard model which can describe the decoherence with 
population decay is given by T^sB = W± with 

n± = Y,f^ {9k±(^- gk±a+ ® 6fc) . (4) 

k 

It has been proposed that decoherence as well as population decay can be suppressed by ap- 
plying periodic tt pulses about the z-axis if the interval At is much shorter than the correlation 
time.'^'^'^ The present authors^^ confirmed this by numerically calculating the time evolution of 
the reduced density matrix for a qubit based on the time-convolutionless projection operator 
approach. Protopopescu^^ reported a numerical study on the subject similar to this. However, 
they employed the approximation in which the commutator ['Hp{t),Ti.±] {Ti.p{t): external pulse 
field) is completely neglected. We have no reason to justify the approximation. 
We focus on the case with 

nsB = nz + H±. (5) 

This case has not been studied extensively. This interaction Hamiltonian TCsb induces not 
only the pure dephasing but also the decoherence with popTilation decay, and describes more 
general situations than those treated in rcfs. 6, 7 and 14. Each ir pulse about the .x-axis (z- 
axis) induces a bit- flip (phase-flip). Periodic bit-flips cancel out the influence of TC^, while 
periodic phase-flips cancel ovit the influence of TC±. Let and Tc± be the correlation times 
for Tiz and 7i±, respectively. It has been shown that if At is much shorter than Tqz and Tc±, 
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the dccohercncc caused by TCsb = T~(-z + T~i-± can be suppressed by a sequence of both bit-flips 
and phase-flips.^' "'^'^ The simplest sequence is composed by the alternate application of a bit- 
flip and a phase-flip. We call it bp sequence. Although its effectiveness has been proved in 
the short interval limit of At ^ t^z and t^-j-, a realization of this condition is not necessarily 
possible in actual experiments. It is of interest to study the effectiveness in experimentally 
relevant cases of At < Tcz and Tc±. Furthermore, we are also interested in whether the bp 
sequence is always more effective than the sequence consisting of only bit-flips. We call the 
latter one bb sequence hereafter. If Tcz ^ Tci, we expect that the bp sequence is much more 
effective than the bb sequence. However, in the regime of t^z <C 7^;^, the bb sequence may 
predominate the bp sequence. 

In this paper, we study the time evolution of a qubit coupled with a quantum environment 
under periodic vr pulses, to examine the effectiveness of the bp and bb sequences when the 
interaction Hamiltonian is given by eq. (5). It is assumed that qubit-environment coupling is 
week. We employ the time-convolutionless projection operator approach following ref. 14, and 
derive the equation of motion which describes the time evolution of the reduced density matrix 
for a qubit. Our approach is simpler than that presented in ref. 12, and more accurate than that 
presented in ref. 13. Besides the approaches in refs. 12 and 13, several approximate methods 
based on the influence-functional formalism may be applicable. However, such methods 
are too complicated for our restricted purpose. The resulting equation has an advantage 
that it is applicable to arbitrary pulse sequences consisting of bit- and/or phase flips. We 
examine the three cases of Tc±/tcz = 2, 5 and 50 for a fixed value of Tc, where Tc is defined by 
Tc"^ = Tcz~'^ + Tc±^'^ . The pulse interval is chosen as At/Tc = 2^'^, 2^'^ and 2^'^. We show that 
the bp sequence is much more efi'cctive than the bb sequence when Tc±/tcz = 2. However, we 
also show that their relative effectiveness can be reversed in the cases of Tc±/tcz = 5 and 50 
depending on At/Tc- Only in the short-interval limit, the bp sequence is always more effective 
than, or at least as effective as, the bb sequence. We conclude that the effectiveness of the 
pulse sequences depends on a relative strength of the two decoherence processes when a pulse 
interval is slightly shorter than Tq. 

2. Formulation 

We derive the equation of motion for the reduced density matrix for a qubit. Let us 
consider the case where the qubit-environment interaction is given by eq. (5) and periodic 
TT pulses with an interval At are applied. We also assume that the duration of each pulse is 
infinitely short. The bb sequence consisting of only bit-flips is described by 



oo 




(6) 
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The bp sequence consisting of both bit- and phase-flips is described by 



oo 

vr 



^bp(i) = 2 Yl ["^(^ ~ ~ cos wo* + <Ty sin wo*) + S^t - 2jAt)a^] (g) Ig, (7) 

i=i 

where a bit-flip is applied at t = (2j — l)At while a phase-flip is applied at t = 2jAt. The total 
Hamiltonian is Ti.{t) = Tio + Hp{t), where 'Hp{t) represents either Hhh{t) or Hbp(i)- Although 
we treat only the bb and bp sequences, the following argument can be applied to any pulse 
sequences consisting of bit- and/or phase-flips. Let p{t) be the density matrix of the whole 
system. It obeys the Liouville-von Neumann equation 

^p{t) = L{t)p{t) (8) 

with 

m- =-^[H{t),.]. (9) 

We decompose H{t) as H{t) = Hi{t) + Hsb, where Hi{t) = Hs+HB + fi-pit). It is convenient 
to introduce p{t) which is defined as 

p(t)= Tae-i'o'^-^i(-)p(t), (10) 

where Tg, is the anti-time-ordering operator and Li{t)- = —^iHiit), •]. Prom eq. (8), we obtain 

^p(t) = LsB(i)p(t), (11) 
where LsB(i)- = -■^[^sB(i), •]• Here, fisBit) is given by 

= raeS^*'^^^p('")7i:sB(i)re-s^*'^^^p(^\ (12) 
where T is the time-ordering operator, np{t) = ct(Ws+WB)t^p(-^)g-i(Ws+-HB)t ^nd 
HsB{t) = e^(^s+WB)t7^gge-i(^s+HB)« 

k 

{9l±<'- ® ble'(-'=-°)* + 9k±a+ bfee-'^-'^-")*) . (13) 

k 

We can show that the time-ordered exponential factor in eq. (12) is obtained as 

Is (8) 1b for < t < At, 

io-a; (g) 1b for At <t < 2At, 

{iax X iax) 1b for 2At <t < SAt, (14) 

{iux X iax X iax) 1b for 3At <t < 4At, 
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Is 1b 
iax ® 1b 
(i(T2 X ia^ 



for < t < At, 
for At <t < 2At, 
for 2At<t< 3At, 



(15) 



(icTa; X icr^ X i(7x) (g) 1b for 3At <t < AAt, 



for ?t:p(t) = Hhpit). Let A(,:(t) (A^2(t)) be the number of bit-flips (phase-flips) within [0,t]. If 
the bb sequence is appHed (i.e., 'Hp{t) = Hhhit)), we find that 

t 



N,{t) = 



At 



0, 



(16) 
(17) 



where [a] denotes the integral part of the real number a. In the case where the bp sequence is 
applied, we find that 

't + Af 



A^x(i) 
N,{t) 



2At 

t 



2At 



If we define 



we obtain 



<Tb(i) 



(7+ if Nx (t) is even, 
(j_ if Nx (t) is odd, 



where 



and 



(18) 
(19) 

(20) 

(21) 
(22) 

(23) 



Equation (21) holds for any pulse sequences consisting of bit- and/or phase-flips, although we 
have assumed that 'Hp{t) is equal to Ti-hhit) or '^hpit)- Note that each bit-flip vr pulse changes 
the sign of Ti.z{t) and each phase- flip vr pulse changes the sign of 'H±{t). The periodic sign 
change due to bit-flips (phase-flips) effectively cancels out the qubit-environment interaction 
Hzit) (H±{t)). This is the reason why the decoherence is suppressed by the vr pulses. 

Our interest is focused on the time evolution of the qubit, which is described by the 
reduced density matrix defined by 

yPoi{t) Poo{t)J ' 



psit) = TTB{m} 



(24) 
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where jjn (poo) represents the population in the upper (lower) state and the off-diagonal terms 
characterize the coherence of the qubit. We derive the equation of motion for psit) following 
ref. 14. In doing so, we assume that the qubit and the environment are uncorrelated at t = 0, 
i.e., 

m = PS(0) ® pb(0), (25) 
and the environment is initially in thermal equilibrium at temperature T, i.e., 

Pb(0) = pb = n (l - 6-'^'^'=) e-^'^'^^l^'' (26) 

k 

with P = 1/kBT {kB - the Boltzmann constant). We introduce the projection operator P which 
is defined as P- = Ttb {•} <8) Pb- We observe that Pp{t) = ps{t) <8) Pb- Adapting the time- 
convolutionless projection operator formalism^^ with P to eq. (11), we derive the equation of 
motion for ps{t), 

^ps(t) = Tre 1^ dhLsB{t)LsB{ti) [ps{t) ® Pb] | , (27) 

where we have neglected higher-order terms by assuming that the qubit-environment interac- 
tion is weak enough. In calculating the right-hand side of eq. (27), we note that the dccohcrence 
without population decay is mainly caused by the modes with u;^ ~ in Tizit) and the deco- 
herence with population decay is mainly caused by the modes with uj^ ~ (jJq in 'H±{t). Thus, 
we approximately neglect the cross terms between Tizit) and 'H±{t) in eq. (27). 
Using eqs. (13) and (21), we obtain 

^Pii(t) = -7ii(i)Pii(i)+mi(i), (28) 
^Pio(t) = -7fo(i)Re[pio(i)]-i7?oWlm[pio(t)], (29) 



where 



711 (i) 



= 2^b.±Pcoth (||,) J^dU 



1 + (_i)iV4t)+JV4ti) 
X — ^— cos{iOk-u}o){ti-t), (30) 

1 + (_l)JVx(t)-HiV4*i) 
X ^ — ^— cos{ujk-u)o)iti-t), (31) 

7fo(i) = 4j^b,.pcoth(|^)jJ|*dti(-l)^^W+^^(*^)cosa;,(ti-t) 

+ El^-^oth(^)^^.M-i)-^*)-^^*^) 

^(•_l)Ar^{t)+Ar^(ti)gi(-i)^(*)K-'^o)(ti-t)^ (32) 
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Fig. 1. Time evolution of pn as a function of t/r^ at k^T/tvjJc = 10~^. The ratio of the qubit- 
environment correlation times is Tc±/tcz = 2. The pulse interval is (a) At/Tc = and (b) 
2~^. The solid lines and dotted lines correspond to the cases of the bb sequence and the bp 
sequence, respectively. The dashed lines represent the case where n pulses are absent. 



In deriving eqs. (28) and (29), we have used the two relations pii{t) + poo{t) = 1 and pio{t) = 
PQi{t). Equations (28) and (29) describe the time evolution of the qubit under tt pulses. It 
should be emphasized that the above equation can be applied to any pulse sequences consisting 
of bit- and/or phase-flips. Introducing the spectral functions Iz{to) and I±{io) defined by 

IzH = Y^\9kzfS{u-uJk), (34) 

k 

4('^) = ^\9k±fS{uj - Uk), (35) 
k 

we rewrite eq. (33) as 

7?o(i) = 4^" du^Uu) coth 1^ dh (-l)^^W+^^(*0 cosa;(ti - t) 

+ rdioUu) coth (^) /*dti(-i)A^^W+^4*i)eK-i)-^«(--o)(ti-t).(36) 

Jo y^KBTj Jo 

Equations (30), (31) and (32) are rewritten in the manner similar to eq. (36). 

3. Numerical Results 

On the basis of the resulting equation of motion, we study the time evolution of the 
reduced density matrix psit) in the presence of periodic vr pulses with an interval At. 
We numerically calculate pn{t), \pio{t)\ and A6{t), where A9{t) is defined by A9{t) = 
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arg{pio(t)} — arg{pio(0)}. We assume that the spectral function is given by 



Iz{uj) = G^ue 



-Ul/U)c 



(37) 



I±{uj) = G±Loe 



(38) 



where cOc is a cut-off frequency, and and G±- are coupHng constants. The following parame- 
ters are employed: ooq/lOq = 0.1 and U-qT/ (fujc) = 0.001. In terms of the coupling constants, the 
correlation times for Tiz and 7i± are given by Tcz = {1/2Gz)^^'^uJc~^ and Tc± = {2/G±)^^'^uJc^^ , 
respectively. The initial behavior of pio{t) is expressed as pio{t) = pio(O) exp[— (t/xc)^], where 



We fix Tc at Tc = 0.4 x (2'k/uJc) and examine the cases of Tc±/tcz = 2, 5 and 50. While the 
pure dephasing is slightly stronger than the decoherence with population decay in the case 
of Tc±Itcz = 2, the former predominates the latter in the case of Tc±/tcz = 50. For each case, 
we calculate the time evolution with At/Tc = 2~^, and 2~^. We choose (|0) -|- |l))/\/2 
and (|0) -I- i|l))/\/2 as initial states. From numerical results, we find that the decoherence 
suppression is relatively weak in the latter case. Thus, we present only the results for (|0) -|- 



In Fig. 1, we display pii{t) with (a) At/Tc = 2^'^ and (b) 2~^ in the case of Tc±/tcz = 2. 
We do not display the results for Tc±/tcz = 5 and 50 because variations of pii{f^ are very 
small. The decrease of pu represents the population decay. We observe that the population 
decay is suppressed by either of the bb or bp sequences. It should be emphasized again that 
the bb sequence cannot suppress the influence of 7i±, which induces population decay, while 
the influence is suppressed by the bp sequence containing phase-flip vr pulses. The suppression 
of the population decay under the bb sequence is simply attributed to periodic exchanges of 
the upper-state population with the lower-state population due to bit-flip vr pulses. 

We compare pio{t) under the bp sequence with that under the bb sequence. In Fig. 2, we 
display |/5io(i)| with At/Tc = 2^^, 2~'^ and 2^^ when (a) Tc±/tcz = 2, (b) 5 and (c) 50. The 
decrease of |/9io(t)| represents the decoherence. Figure 2 shows that in the case of Tc±/tcz = 2, 
the decoherence suppression under the bp sequence is much more effective than that under 
the bb sequence. For the bb sequence, the decrease of At from At/Tc = 2~^ to 2^^ does not 
result in a notable improvement. This fact is attributed to a shortcoming of the bb sequence 
which cannot cancel out the influence of TC±. However, the effectiveness of the bb sequence is 
improved with increasing Tc±/tcz since the influence of Ti^ is relatively reduced. In the case 
of Tc±/tcz = 5, the bb sequence is more effective than the bp sequence for At/Tc = 2~^, and 
is predominated by the bp sequence for At/Tc = 2~^ and 2~^. In the case of Tc±/tcz = 50, 



i|l))/V2 with 




(40) 
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Fig. 2. Time evolution of \pio\ in the case of (a) Tc±/tcz = 2, (b) 5 and (c) 50. The sohd hncs and 
dotted Unes correspond to the cases of the bb sequence and the bp sequence, respectively, with 
^t/Tc = 2~^, 2~^ and 2~^ from lower to upper lines. The dashed lines represent the case where tt 
pulses are absent. 



the bb sequence is more effective than the bp sequence for At/rc = and 2~^, and is as 
effective as the bp sequence in the shortest interval case of At/Tc = 2~^. The above results 
indicate that the effectiveness of the pulse sequences depends on a relative strength of the two 
decoherence processes especially when a pulse interval is slightly shorter than Tq. Only in the 
short-interval limit, the bp sequence is more effective than, or at least as effective as, the bb 
sequence regardless of Tc±/tcz- 

The above results are explained as follows. The bp sequence contains both phase-flips and 
bit-flips, and thus can cancel out the influence of both Hz and H±, while the bb sequence can 
cancel out the influence of only Hz- Thus, the bp sequence is much more effective than the 
bb sequence in the case of Tc±/tcz ~ 1, where the influence of both Hz and H± is important. 
However, the interval of two adjacent bit-flips in the bp sequence is equal to 2At, and is twice 
longer than that in the bb sequence. This indicates that the bp sequence reduces the influence 
of Hz less effectively than the bb sequence. Thus, in the case of S> 1, where Hz is more 

important than H± , the bp sequence becomes less effective than the bb sequence unless At is 
very much shorter than t^z- 

In Fig. 3, we display the phase shift AO{t) for At/rc = and in the case of 
Tc±/tcz = 2. Obviously, AO{t) = in the absence of qubit-environment interactions. We 
observe that the phase shift is suppressed by either of the bb or bp sequence. The phase shift 
in the cases of T^^i/ij^ =5 and 50 is very small, so we do not display the corresponding results. 

4. Summary and Conclusion 

To study the effectiveness of the pulse control, we have studied the time evolution of a 
qubit coupled with a quantum environment under a periodic vr pulses. We have considered 
a general interaction Hamiltonian which induces not only the pure dephasing but also the 
decoherence with population decay. We have examined the effectiveness of the bb sequence 
consisting of only bit-flips and the bp sequence consisting of bit- and phase-flips. We have 
derived the equation of motion for the reduced density matrix. The resulting equation is 
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A^(^)[rad] A^(f)[rad] 




0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 

(a) t/rc (b) t/rc 

Fig. 3. Phase shift Ae{t) in the case of %±/%z = 2 with (a) At/rc = 2'^ and (b) 2"^. The solid Unes 
and dotted hnes correspond to the cases of the bb sequence and the bp sequence, respectively. 
The dashed lines represent the case where tt pulses are absent. 

applicable to arbitrary sequences consisting of bit- and/or phase-flips. 

By numerically solving the equation of motion, we have shown that the effectiveness of 
the pulse sequences depends on a relative strength of the two dccohcrence processes especially 
when a pulse interval is slightly shorter than t^. We have found that the bp sequence is 
more effective than, or at least as effective as, the bb sequence regardless of Tc±/tcz only in 
the short-interval limit. This means that if we cannot prepare periodic vr pulse with a short 
enough interval, the bp sequence is not necessarily the best choice unless Tc±/tcz ^ 1. 
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